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Abstract 

We present non-local dynamics of Bell states in separate cavities. It is demonstrated that (i) 
the entanglement damping speed will saturate when the cavity leakage rate 7 > 0.4; (ii) the 
synchronism relationship between the fidelity and the concurrence depends on the initial state; 
(iii) if the initial state is l/v2(|01) + 1 10) ) , the dynamics of entropy is opposite to that of fidelity. 
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I. INTRODUCTION 

In contrast with the extensively investigated static entanglement 
dynamic entanglement under the influence of variant environments is one of the most 



important and largely unexplored problems in the fie 



d of quantum teleportation, quantum 



computation and quantum communication [3, [s[ 0, 1(3]. It is not only involved with the 
foundation of quantum mechanics, but also a fundamental issue in creating, quantifying, 
controlling, distributing and manipulatin g th e entangled quantum bits, which are composed 
of spin-1/2 atoms in differe nt p roblems [ll|, ll2j, |l3j . An entangled system is in such a state 



that cannot be factorized 



14j in its Hilbert space. And the most familiar and widely used 



examples are Bell states. The two particles or atoms of spin-1/2 are correlated no matter 
how long distance is between them. Generally, due to different kinds of quantum reservoir, 
the entanglement degree between them vanishes asymptotically. However, if the reservoir 
consists of, e.g., only one or two electromagnetic field modes, then the entanglement may 
decrease abruptly and non-smoothly to zero in a finite time [9, [jj], [if]], which is a new 
nonlocal decoherence called entanglement sudden death (ESD). Therefore, demonstration 
of the dynamics of Bell states [l^] would have profound implications for understanding of 
the physics in the realization of qubits in experiments. 



So far in quantum optics experiments, Bell states can be generated with trapped ions 



[l8l ] and in cavity quantum electrodynamics (CQED) [ji], |20j], whirl) 
attention. Based on cavity QED systems, schemes (to see Refs. [21 



las attracted much 



22 



23 



24 



25j and 



references therein) have been proposed to implement quantum communications or engineer 
entanglement between atoms in distant optical cavities. In most of them, two separated 
cavities are connected via some channels, for instance, an optical fiber 25]. And in a 



recent paper, Yin and Li 



261 ] investigated a system consisting of two single-mode cavities 



connected by an optical fiber and multiple two-level atoms trapped in the cavities. They 
show that ideal entangling can be deterministically realized between the distant cavities. 
Besides, utilizing a system of two-atoms and two photon modes, Masood and Miller 27] 



used the Jaynes and Cummings model 



28 



appropriate models for exploiting the dynamics of entanglement 



which is considered to be one of the most 

3211. in the 



29 



30, 



31 



rotating wave approximation to study entanglement of more than one atom with vacuum. 
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The photon modes in their model are uncoupled, however, the leakage of cavities and the 
effect of temperature (Yet recently, in other models the two-qubit entanglement dynamics 
for a finite-temperature environment has been discussed in 0, 34 1) are actually not 
considered thoroughly. 



In this paper, we consider a quantum model with two identical two-level atoms or pseudo- 
spins of 1/2 (as an open subsystem with qubits labelled si and s^) and two single-mode 
cavities (labelled 1 and 2 correspondingly). The atom Sj (j = 1 or 2) is embedded in and 
coupled only with the cavity mode j, which could be regarded as its bath or environment. 
The two cavities are so far departed that there is no direct interaction between them as well 
as the two atoms. Initially, the two qubits are prepared as a most-entangled states (Bell 



states). The focus o: 
the concurrence 35 



interest is their degrading quantum evolution, which are measured by 



361 ] , the fidelity 37| and the entropy exchange 



38 



391 ]. The calculations 



and physical arguments will be carried out in two conditions: (i) there is leakage of photons 
for the cavities, which are in the vacuum states from the beginning; (ii) the cavities are so 
perfect that the loss of photons from them could be neglected and the two single modes are 
initialled in a thermal equilibrium state with the same temperature. The rest of this paper 
is organized as following. In Sec. [TT] we begin with the model Hamiltonian and its analysis 
derivation; and then we introduce the numerical calculation procedure about the evolution 
of the reduced matrix for the subsystem. Detailed results and discussions can be found in 
Sec. IIHI We will make a conclusion in Sec. [TV] 



II. MODEL AND METHOD 

The master equation for a two-level atom in a single-mode cavity [40|], as one of the two 
partitions in our model, can be taken as 

For density matrix pj, j refers to Si or s 2 ; for the mode operator dj or at, j (1 or 2) represents 
the photon mode coupling with the corresponding atom, jj is the leakage rate of photons 
from the cavity j. Hj describes the Hamiltonian for a subsystem of one atom and one cavity 
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(J = 1,2): 



Hj = -j-a) + (1 + <jUj<>'<> i + gjUjia] + a,)a':. (2) 



where ujj is the energy level difference of atom Sj in cavity j. ej is the detuning parame- 
ter measuring the deviation of the photon j energy from ujj. gj is introduced as another 
dimensionless parameter which suggests the coupling strength between qubit Sj and mode 
j. The x and z components of o are the well-known Pauli operator. The two qubits are 
embedded in remote cavities without direct interaction. Therefore the whole Hamiltonian 
for this two-atom-two-cavity problem is 

H = H 1 + H 2 . (3) 

The whole state of the total system is assumed to be separable before t = 0, i.e. 

p(0) = Ps (0)®p b (0), (4) 

Ps (o) = (5) 

Pb(0) = Pw(0)®p w (0). (6) 

The initial state |^(0)) for the two qubits is one of the Bell states. And the two cavities 
are in their (i) vacuum states Pfej(0) = |0j)(0j| (in this case, we will consider jj ^ 0) or 
(ii) thermal equilibrium states /%(0) = e~ HB ^ kBT / Z (in this one, we set 7 to be zero to 
distinguish the effect of temperature from that of 7), where Hb is the pure bath part of 
the whole Hamiltonian and Z = Tr (e~ HB / kBT ^j is the partition function and the Boltzmann 
constant ks will be set to 1 for the sake of simplicity. 

For the former case, Eq. [T]will be exploited to calculate p{t). For the latter one, Eq. [1] 
is reduced to 

pit) = exp(-iHt)p(0) exp(iHt). (7) 

To determine the dynamics of the density matrix for the whole system, two factors need 
to be considered. The first one is the expression of the thermal bath state. In numerical 
calculations [4l|, we have to expand Pbj{0) {j = 1, 2) to a summation of its eigenvectors with 
corresponding weights determined by its eigenvalues: 

Eg— E m j/T 
\<j>mj)u mj ((j) m j\, UJ mj = (8) 

m. J 



Then for the two single-modes, we have 

Pfel(O) ®Pfe 2 (0) = ^ |0 m l)|0n2)^mn(0n2|(0ml|, ^r, 



-(E ml +E n2 )/T 
Z1Z2 



(9) 



where the subscripts m and n refer to mode 1 and 2 respectively. The second important 
factor is the evaluation of the evolution operator U(t) = exp(iHt). A polynomial expansion 



scheme proposed by us in Ref. 

U(t) = 



42 



43 



441 ] is applied into the computation, 

o+l 00 / .. \ k 



1 + it 



k=0 



it 



+ it 



(10) 



L^(H) is one type of Laguerre polynomials as a function of H, where a (—1 < a < 00) 



distinguishes different types of the Laguerre polynomials and k is 
scheme is of an efficient numerical algorithm motivated by Ref. 



;he order of it. The 



45 



461 ]. which is pretty 



well suited to many quantum problems, open or closed. Additionally, it could give results 
in a much shorter time compared with the traditional methods under the same numerical 
accuracy requirement, such as the well-known 4-order Runge-Kutta algorithm. After the 
density matrix p(t) for the whole system is obtained, the reduced density matrix psif) for 
the two atoms can be derived by tracing out the degrees of freedom of the two single-mode 
cavities. 



III. SIMULATION RESULTS AND DISCUSSIONS 

We discuss three important physical quantities which indicate the time evolution of the 
subsystem, (i) The concurrence. It is a very good measurement for the intra-entanglement 
between two qubits and monotone to the quantum entropy of the subsystem when the 
subsystem is in a pure state. It is defined as: 

C = max{Ai - A 2 - A 3 - A 4 , 0}, (11) 

where A, are the square roots of the eigenvalues of the product matrix ps(cr y ®cr y ) p* s (a y ®cr y ) 
in decreasing order, (ii) The fidelity. It is defined as 

F(t) = Tv s \p ideal (t)ps(t)}. (12) 

where pideaiit) represents the pure state evolution of the subsystem only under Hs, without 
interaction with the environment. In this study, Hs = + ^o^. The fidelity is a 

5 



measurement for decoherence and depends on pideai- It achieves its maximum value 1 only if 
Ps(t) equals to Pideai(t). (iii) The entropy exchange En is defined as En = — Tr(p 5 log 2 ps)- 
It is the von Neumann entropy of the joint state of the subsystem as composed of the two 
qubits in our model. It measures the amount of the quantum information exchange between 
the subsystem and the environment. For the subsystem consisted by two two-level atoms 
(its Hilbert space is 4 x 4), the entropy maximum is log 2 (4) = 2.0. When it reaches its 
maximum value, it means all the quantum information is cast out of the subsystem or the 
quantum subsystem degenerates to a classical state. 

A. Dynamics at different 7 




5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 

CO t CO t 



(a)C(i) (b)-Fd(t) 

FIG. 1: Time evolution for (a) Concurrence, (b) Fidelity with the subsystem starting from 
l/\/2(|00) + 1 1 1) ) at different values of anisotropic parameter: 7 = (solid curve), 7 = 0.2 (dot 
dashed curve), 7 = 0.4 (dashed curve), 7 = 0.8 (dotted curve). The two cavities are initialed as 
|0)i|0> 2 . 

In order to discuss the effect of 7 (We suppose the two cavities have the same loss degree: 
7i = 72 = 7-) and T, all the other parameters are fixed for the sake of simplicity and 
without loss of generality: 

<jj\ = 002 = <jj = 0.4, 
ei = £2 = e = —0.5, 
9i = 92 = 9 = 0.2. 
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5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 
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(a)C(t) (b)Fd(t) 

FIG. 2: Time evolution for (a) Concurrence, (b) Fidelity with the subsystem starting from 
l/\/2(|01) + 1 10)) at different values of anisotropic parameter: 7 = (solid curve), 7 = 0.2 (dot 
dashed curve), 7 = 0.4 (dashed curve), 7 = 0.8 (dotted curve). The two cavities are initialed as 
|0)i|0> 2 . 
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(a)l/V2(|01) + |10)) (b)l/V2(|00) + |11)) 

FIG. 3: Time evolution for entropy of the subsystem from (a) 1/a/2(|01) + |10)), (b) 1/^/2(1 00) + 
1 11)) at different values of anisotropic parameter: 7 = (solid curve), 7 = 0.2 (dot dashed curve), 
7 = 0.4 (dashed curve), 7 = 0.8 (dotted curve). The two cavities are initialed as 1 0) 1 1 0)2 - 

And we choose l/-\/2(|00) + 1 1 1 > ) and l/v^(|01) + |10)) as two different initial states for 
the subsystem. 



We first discuss the effect of the photon loss rate 7. It is evident that with a larger 7, 
the entanglement degree of the subsystem will decrease in a faster speed, which could be 



verified by Fig. 1(a) and Fig. 2(a). The tendency of the two cases is similar. If 7 = 0, 



the concurrence will oscillate periodically with time and will not be dissipated; but the 
peak value of it will never reach 1.0. On the whole, the curves of the concurrence are not 
perfectly harmonic, which is a little different with the results gotten in previous works. It is 
due to the stochastic and irrelevant microscopical processes (the spins drop from the excited 
state by emitting a photon or jump to the excited state by absorbing a phonon) inside the 
two different cavities. And the dynamics of the concurrence stems from such numerous 
processes, so the evolution is approximately harmonic but not perfect. When 7 > 0, the 
concurrence drops abruptly to zero in a short time. It coincides with the description about 



entanglement sudden death (ESD) in Ref. [l5|, |47| that "after the concurrence goes abruptly 
to zero, it arises more or less from nowhere". This is an example of ESD. The photons 
leaking out of the cavities greatly reduced the nonlocal connection between the two qubits. 
When 7 is bigger than 0.4, the speed of ESD is saturated and we almost cannot distinguish 
the curve of 7 = 0.4 from that of 7 = 0.8. The state l/v^QOl) + |10)) seems more robust 
than l/v2(|00) + |11)). It is verified that, for example, in the condition of 7 = 0.2, the 



concurrence of the former state (to see the dot dashed line in Fig. 2(a)) decreases to zero 
at ut = 17.472 for the first time, while the latter one does at ut = 15.168. 



The fidelity dynamics of the two Bell states, however, is much different from each 



other. In Fig. 1(b), when there is no cavity leakage, the fidelity is not synchronous with 



the concurrence. Concretely, along that curve, only when ut = 37.376, most component 
of the subsystem state is recovered to its initial one with high fidelity F(t) = 0.944346. 
While the other three peaks inside the interval of ut G [5.0, 35.0] along the solid line 



in Fig. 1(a) are just fake phenomena: although the entanglement degree between the 



two subsystem atoms is high, but the state of the subsystem is different from the initial 
one. The tendency of 7 = 0.2 and 7 = 0.4 is in a similar manner: the curves decrease 
smoothly and move towards F = 0.3 ~ 0.4. While the curve of 7 = 0.8 has an oscillating 
dynamics with gradual shrinking amplitude. Finally, its value approaches to F = 0.35. 
In Fig. 2(b), however, when 7 = 0, the fidelity is synchronous with the concurrence. 



Thus the fidelity of the state 1 01) + 1 10) is more robust than that of 1 00) + 1 1 1 ) under 
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the same environment. These contrasts between the two initial state can be noticed 
in other works. It implies the physics essence of them is different, although both of 
them are of the most entangled states. Yet the other three cases with leakage 7 > 
have almost the same kind of dynamics for the two states. The descend speed of the fi- 
delity decreases with time and the value of fidelity approaches to about F = 0.2 at ut = 40.0. 



The quantum evolution of entropy exchange En is depicted in Fig. [31 Initially, the entropy 
of the two entangled qubits equals to 0, which means that all the quantum information is 
kept in the entanglement between the two qubits. We find that when 7 = 0, the evolutions 
of the two Bell states are a little different but the oscillation periods of them are almost the 
same. Although there is no leakage, but the bath, the two single-mode cavities, will absorb 
some of the quantum information inside the subsystem, which means the entropy never goes 



back to zero as initialled. To compare the solid line in Fig. 3(b) (Fig. 3(a)) with that in 



Fig. 1(a) (Fig. 2(a)), we notice that the tendency of the entropy is opposite to that of the 
concurrence. The dynamics difference from the two initial states can almost be removed by 



introducing non-zero 7 as it is shown in the comparison of Fig. 3(a) and Fig. 3(b) It is 
evident that with a larger 7, more quantum information of the subsystem is transferred into 
the bath. 



B. Dynamics under different T 

In this subsection, we turn to the effect of temperature of the cavity modes. It is found 
that when the temperature is comparatively low, T = 0.25u, the entanglement degree 



oscillates but will not corrupt to a sudden death (to see the dot dashed line in Fig. 4(a) and 



Fig. 5(a)). When it goes up to a moderate temperature T = 0.5a;, the ESD happens. The 
first moment at which the concurrence decreases to zero and does not revive immediately 
is dependent on the initial state. For l/v2(|00) + |11)), it takes place at cut = 4.192, which 
is much earlier than that of l/\/2(|01) + |10)), cut = 11.104. This difference coincides with 
the comparison in subsection IIII Al which means the entanglement damping speed of the 
state |01) + 1 10} is slower than that of state 1 00) + |11)). It is hinted that the bath influence 
on the former state is comparatively inapparent. Yet both of them can still revive to a 
certain extent 0.2 ~ 0.3 after some time. When the temperature increases to T > 0.75u;, 
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5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 

CO t CO t 



(a)C(t) (b)Fd(t) 

FIG. 4: Time evolution for (a) Concurrence, (b) Fidelity with the subsystem starting from 
l/\/2(|00) + |11)). The two cavities are initialed in thermal states at different temperature: T = 
(solid curve), T = 0.25a; (dot dashed curve), T = 0.5a; (dashed curve), T = 1.0a; (dotted curve). 




5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 

CO t CO t 



(a)C(t) (b)-Fd(t) 

FIG. 5: Time evolution for (a) Concurrence, (b) Fidelity with the subsystem starting from 
l/\/2(|01) + 1 10)). The two cavities are initialed in thermal states at different temperature: T = 
(solid curve), T = 0.25a; (dot dashed curve), T = 0.5a; (dashed curve), T = 1.0a; (dotted curve). 

it can not revive in the future after the first sudden death happens for both Bell states. 
Till the temperature is as high as T = 1.0a;, the concurrence of both cases falls with a very 
quick speed to zero. Obviously, the temperature will destroy the initial most-entangled 
states even if it is much lower than the energy bias ui of the two-level atoms. After ESD 
happened, the quantum oscillation from the local thermal bath may help to entangle the 
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cot cot 



(a)l/V2(|01) + |10)) (b)l/V2(|00) + |11)) 

FIG. 6: Time evolution for entropy of the subsystem from (a) l/\/2(|01) + |10)), (b) 1/V2Q00) + 
1 11)) on the condition that the two cavities are initialed in thermal states at different values of 
temperature: T = (solid curve), T = 0.25a; (dot dashed curve), T = 0.5a; (dashed curve), 
T = 1.0a; (dotted curve). 

two qubits, but this positive effect is neglectable when the temperature is high enough. 
Then the entanglement between the Bell states is damped forever. 



Under variant temperature, the different dynamics of fidelity, which rely on the initial 
state, are shown in Fig. 4(b) and Fig. 5(b). In a short interval after ut = 0, higher 
temperature means faster damp speed for both initial states. Yet in a long time scale, their 
actions are totally different. For T < 0.5a;, the four curves evolves pseudo-periodically, 
but the period of l/>/2(|00) + |11)) is much larger than that of 1/V2(|01) + |01)). In Fig. 



4(b), the dotted curve of T = 1.0a; fluctuates with time, whose amplitude damps from the 



beginning time and gets some revival when ut > 25.0. In Fig. 5(b) , the oscillation evolution 
at temperature T = 1.0a; is in the same manner as those at T = 0.25a; and T = 0.5a;. 
And their amplitudes and peak values of the fidelity decrease with increasing temperature. 
Obviously, an environment at higher temperature destroys the fidelity of the subsystem 
even stronger. 



In Fig. [61 we give the dynamics of the entropy exchange under different temperatures. 
Although there are some differences between the two sub-figures even when T is as high as 
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0.5a;, but the tendency of the two cases are almost the same. With higher T, the entropy 
increases with faster speed and behaves an oscillation evolution with smaller amplitude. It 



is important to find that the entropy exchange in Fig. 6(a) exhibits an opposite behavior 



in comparison with that of the fidelity in Fig. 5(b) The periods of the two evolutions 
are the same and when the fidelity experiences a peak value, (For instance, for the dot 
dashed curves in the two figures, when T = 0.25a;, the four peaks appear at ut = 8.160, 
ut = 18.656, ut = 27.776 and out = 37.696 in the give n interval) the entropy is at the 
corresponding valley point and vice versa. In Ref. 48j, the authors found the entropy 
exchange exhibit the behavior opposite to that of the concurrence. That coincides with 
our results. In our model, when the temperature is not too high and the initial Bell state 
is l/v^dOl) + 1 10) ) , the concurrence and fidelity evolve synchronously, then the dynamics 
of entropy is also opposite to that of the concurrence. However, when the temperature is 
high enough, the ESD will make this synchronization relationship invisible. So the opposite 
relationship between the entropy and concurrence is lost. Therefore we have to conclude 
that the more quantum information of the subsystem transfers to its bath, the more fidelity 
of that decreases during the time interval if \ip(0)) = l/v2(|01) + 1 10) ) . 



IV. CONCLUSION 



In conclusion, we investigate the dynamics of two distinct uncoupled qubits embedded 
respectively in two single-mode cavities with leakage, which constitute the environment in 
our model. The subsystem consisting of the two qubits is initially prepared as one of the 
Bell states and the cavities as vacuum states or thermal equilibrium states. Under these two 
situations, the concurrence, the fidelity and the entropy exchange are used to portrait the 
subsystem dynamics from different views. Polynomial expansion method is applied into the 
numerical calculation. It is found that (i) for the leaky cavities, when the loss rate 7 > 0.4, 
the speed of entanglement sudden death achieves its maximal value; (ii) The evolution of 
different Bell states can easily be distinguished by their fidelity dynamics; (iii) the behavior 
of entropy exchange is opposite to that of the fidelity and concurrence in certain conditions. 
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